Accuracy thresholds of topological color codes on the hexagonal and square-octagonal 

lattices 
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Accuracy thresholds of quantum error correcting codes, which exploit topological properties of 
systems, defined on two different arrangements of qubits are predicted. We study the topological 
color codes on the hexagonal lattice and on the square-octagonal lattice by the use of mapping 
into the spin glass systems. The analysis for the corresponding spin glass systems consists of the 
duality, and the gauge symmetry, which has succeeded in deriving locations of special points, which 
are deeply related with the accuracy thresholds of topological error correcting codes. We predict 
that the accuracy thresholds for the topological color codes would be 1 — p c — 0.1096 — 8 for the 
hexagonal lattice and 1 — p c — 0.1092 — 3 for the square-octagonal lattice, where 1—p denotes the 
error probability on each qubit. Hence both of them are expected to be slightly lower than the 
probability 1 — p c — 0.110028 for the quantum Gilbert- Varshamov bound with a zero encoding rate. 



I. INTRODUCTION 

In quantum information processing, we encode infor- 
mation into a quantum state. However the quantum 
state suffers from decoherence by surroundings. We need 
an encoding technique to restore errors on the quantum 
state by decoherence, namely the error correcting code. 
It has been pointed out that error correction of the quan- 
tum state has a deep connection with spin glass models 
and lattice gauge theories in finite dimensions. An out- 
standing instance has been found between the random 
bond Ising model, as well as the random plaquette gauge 
model, and the topological toric code [l| . Recently it has 
been shown that a more useful quantum error correcting 
code, the topological color code, also has a connection 
with the spin glasses [2|, y] . 

These fascinating bridges between the spin glasses and 
the quantum error correcting codes are established in 
a special subspace in spin glasses. This subspace is 
known as the Nishimori line which goes through the 
phase boundary between the ferromagnetic and param- 
agnetic phases as shown in Fig. [1] 4, |5|. The ferromag- 
netic phase on the Nishimori line represents the region 
where errors on the quantum state can be restored, and 
vice versa. Thus the critical point on this subspace in 
spin glasses corresponds to an accuracy threshold of the 
quantum error correcting code. The critical point on this 
line is called the multicritical point. There are little exact 
and rigorous approaches for spin glasses in finite dimen- 
sions. One of the powerful theoretical techniques is the 
gauge symmetry on the Nishimori line. On this line, we 
can calculate the exact value of the internal energy, eval- 
uate the upper bound for the specific heat, and obtain 
some sets of rigorous and exact inequalities by the gauge 



some sets ot rigorous and exact inequalities by trie gauge 
symmetry 0, E| ■ I n addition, the analytical theory esti- 
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FIG. 1: Schematic picture of the phase diagram of the ±J 
Ising model on two-dimensional lattice (left panel) and on 
higher dimensions (right panel). The vertical axis expresses 
the temperature T, and the horizontal line denotes the con- 
centration p of the antiferromagnetic interactions. The mul- 
ticritical point (MCP) is described by the black point. The 
Nishimori line is drawn by the dashed line. For higher di- 
mensions, not only the ferromagnetic (FM) and paramagnetic 
phases (PM) but also the spin glass phase (SG) exists. The 
phase boundary under the Nishimori line is believed to be 
vertical by the argument of the gauge transformation [J. 



mating the precise location of the multicritical point has 
been established by a technique with the duality and the 
replica method [y, 0, H, HI, LLC} . In the present study, we 
analyze two spin glass models, the random three-body 
Ising model (± J type) on the triangular lattice and the 
Union- Jack lattice. These models have deep connections 
with the topological color codes defined on the hexagonal 
and square-octagonal lattices, which implements several 
useful quantum gates 0, HJ. We present the analysis 
of estimations of the error thresholds of these topolog- 
ical color codes through the analysis of the location of 
the multicritical point for the random three-body Ising 
model on two lattices. As a result, we predict the accu- 
racy threshold of the color code on the hexagonal lattice, 
which corresponds to the random three-body Ising model 
on the triangular lattice, would be 1 — p c = 0.1096 — 8 
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FIG. 2: Triangular lattice and Union-Jack lattice. The site 
variable <j) is located at each site, and the plaquette variable 
ki is introduced at each triangular face. The plaquette vari- 
ables form a structure of the hexagonal and square-octagonal 
lattices. 



and that on the square-octagonal lattice, which has a 
connection with the random three-body Ising model on 
the Union- Jack lattice, would be 1 — p c = 0.1092 — 3, 
where 1—p represents the error probability on qubits by 
decoherence. 



II. DUALITY 

Let us review several known facts in the present section 
to fix the notation and prepare for the developments in 
the next section. We take the three-body Ising model on 
the triangular lattice as an example here for simplicity. It 
is straightforward to generalize the following arguments 
to the Union- Jack lattice. 

The Hamiltonian of the three-body Ising model on the 
triangular lattice is 



H 



-J ^2 ViVjVk, 

(ijk) 



(1) 



where the summation is taken over all the triangles (both 
of up-pointing and down-pointing triangles), and a is the 
Ising spin variable taking ±1 on sites as in Fig. [21 This 
model is known as a self-dual model [ill, E3 , and has a 
single transition point T c = 2.26919 which is the same 
as the ordinary two-body Ising model on the square lat- 
tice [13j. We use another expression for convenience as 
follows, 



H 



J COs{7r( 

(ijk) 



fa +(f>k)}- 



(2) 



Here the spin variable is changed into <j) taking and 1. 
It will be useful to review the duality transformation for 



the non-random three-body Ising model on the triangular 
lattice through the Wu-Wang formalism 14], which is 
more suitable for application to spin glasses. We consider 
the face Boltzmann factor for an elementary triangle, 



A J. , J. | i — A'cOsfTT^i+l^+^fc)} 
■^Ipi+cflj+tpk — C ) 



(3) 



where K is the coupling constant K — (3 J. We introduce 
two-component Fourier transform of the face Boltzmann 
factor as, 



h a/2 V 



e K cos irk) 



(4) 



If we use this quantity and another variable fc/, which 
takes or 1, on each elementary triangle as in Fig. [21 
we can rewrite the face Boltzmann factor as, similarly to 
case for the ordinary two-body Ising model on the square 
lattice 0], 
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1 

V2 



fei=0,l 



(5) 



Using this expression, we can rewrite the partition func- 
tion of the three-body Ising model on the triangular lat- 
tice as, 



1 
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i}{k,} I 



(6) 

i+4>3+<j>k) (y\ 



where N is the number of triangles. We take the sum- 
mation over (j), and obtain constraints expressed by the 
Kronecker delta with modulus 2 on each site i on the tri- 
angular lattice 28(ki + k% + ks + + fes + ke), where ki 
surrounds each site as in Fig. O The partition function 
is then rewritten as, 



Z{K) 



N 



JJ(5(fci + k 2 + k 3 + fc 4 + k 5 + fe 6 ). (8) 



We introduce another variable fc/ = <pi + (fj + tp^(vaod 
2), which satisfies all the constraints, and obtain 



(9) 



where ipi also takes and 1 and is located at each site. We 
obtain two expressions for the partition function Z(K) by 
the original and dual face Boltzmann factors as in Eqs. 
^ and We extract the principal Boltzmann factors 
Ao and Aq as, to measure the energy for both of these 
expressions from the state with all-spin up, 



A* N z(ut), 



(10) 
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where z stands for the normalized partition function 
z(ut) = Z/A* and z{u\) = Z/{A*) N . The quantities u t 
and u* are defined as u\ = A\/Aq and u* = A*/Aq. Each 
partition function is now reduced to a single-variable 



function of 



Ml 



-2 if 



and u\ — tanh-ftT. We can identify 



the critical point as a fixed point of the duality by solv- 
ing e~ 2Kc — tanhA'c, under the assumption of a unique 
transition. We remark that, at this fixed point, an ap- 
pealing equation for two principal Boltzmann factors is 
satisfied, Aq = Aq. 



III. MULTICRITICAL POINT 

We step in the analysis of the location of the multi- 
critical point for the random three-body Ising model on 
the triangular lattice, which is related with the accuracy 
threshold of the topological color code on the hexagonal 
lattice [2|. The Hamiltonian is slightly modified as 



H 



J ^ T vk cos( 

(ijk) 



+ (f>j + <t>k), 



(11) 



where Tijk denotes the quenched random coupling. 
Though various types of distribution for r^k can be con- 
sidered, we here restrict ourselves to the ± J Ising model, 
which has a connection with the topological color code. 
The distribution function for the ± J Ising model is 



P(i~i 3 k) = pS(n jk - l) + (i -p)8(Tijk + 

-p)h 



2 cosh K p 
(12) 

. The Nishimori 
,1- 



where K p is defined by cT 2K v = (1 
line is given by the condition K = K p 

We apply the replica method to the ±J Ising model 
The n-replicated partition function after the configura- 
tional average is 



Zu 



II II exp {JfTij-fcof cos( 

{o-i} {ijk) a=l 



+ fa + 4> k )} 



where n stands for the replica number and the angular 
brackets denote the configurational average. We gener- 
alize the duality argument to the n-replicated ±J Ising 
model [1, 0, H, B El- For this purpose it is useful to 
define the face Boltzmann factor Ak (k — 0, 1, • ■ • ,n), 
which represents the configuration-averaged face Boltz- 
mann factor for interacting spins with k triplets giving 
4>i + 4>j + (j>k = 1 (mod 2) among n triplets for a triangle. 
The duality gives the relationship of the partition func- 
tions with different values of the face Boltzmann factor 
as given by 



Z n (A ,Ai,--- ,A n ) 
= Z n (AQ,Al, • • • , A*), 



(14) 



The dual face Boltzmann factors A* k are defined by the 
discrete multiple Fourier transforms of the original face 



Boltzmann factors, which are simple combinations of plus 
and minus of the original Boltzmann factors as the non- 
random case shown above. Two principal Boltzmann 
factors are important pieces of information and given as 
A Q = 2cosh{(n+ 1)K}, and A* = 2™/ 2 cosh™ K, sim- 
ilarly to the case for the random bond Ising model on 
the square lattice @, . We extract these principal face 
Boltzmann factors from Eq. (|14p to measure the energy 
from the all-parallel spin configuration, 



A N z r . 



(ui,u 2 , • • • ,U n ) 



A * N / * * 



;)• 



(15) 



and z n (u\ 1 



are defined as Z u /Aq 



iV 



where z n (ui, ■ ■ 
and Z n /(A* ) N , 

The duality identifies the critical point under the as- 
sumption of a unique phase transition. It is given as the 
fixed point of the duality and is known to yield the ex- 
act critical point for a simple ferromagnetic system as 
reviewed in the previous section. In order to obtain the 
multicritical point of the present replicated spin glass sys- 
tem, we set K = K pi which defines the Nishimori line on 
which the multicritical point is expected to lie. Since 
z n is a multivariable function, there is no fixed point 
of the duality relation in the strict sense which satisfies 
n conditions simultaneously, u\(K) — u\{K) 1 u 2 {K) — 
i*2 (K), • • • , u n (K) — u^(K). This is in sharp contrast to 
the non-random case, in which the duality is a relation 
between single- variable functions. We nevertheless set a 
hypothesis that a single equation A — Aq gives the lo- 
cation of the multicritical point for any replica number 
n [Si LZ|j similarly to the non-random three-body Ising 
model at the critical point. The quenched limit n — > 
for the equation Aq = Aq yields [1, 0] , 

-plogp- (1 - p) log(l - p) = ^log2. (16) 

The solution to this equation is p c = 0.889972. This 
estimation p c — 0.889972 is equivalent to the probability 
for the quantum Gilbert- Varshamov bound with zero rate 
encoding l-p c = 0.110028 0,(3, and consistent with 
the very recent numerical estimation p c — 0.891(2) by 
Katzgraber et al. 0. 

For the random three-body Ising model on the Union- 
Jack lattice, the same duality analysis can be applied and 
we obtain the same solution. In the context of the quan- 
tum error correction, these results mean that the accu- 
racy thresholds for the color codes on the hexagonal and 
square-octagonal lattices are given by 1 — p c = 0.110028 
for error probability on qubits. This consequence would 
not be correct, since two color codes have different com- 
putational capabilities [l7| . For the topological color 
code on the square-octagonal lattice, it is possible to im- 
plement the whole Clifford group of unitary gates gen- 
erated by the Hadamard gate, the {7r/8} gate, and the 
controllcd-NOT gate; while, for that on the hexagonal 
lattice, the {7r/8} gate cannot be implemented. In spin 
glasses, the quenched randomness yields frustration de- 
pending on geometry of lattices, and the locations of the 
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clusters must be chosen to cover the whole lattice under 
consideration as in Figs. [3J and 21 where six examples 
for the triangular lattice and three ones for the Union- 
Jack lattice of the choices of clusters are depicted. It 
is expected that, if we deal with clusters of larger sizes, 
the improved method shows systematic improvements to- 
ward the exact answer on the location of the multicritical 
point P3. 

Then the duality relation for the replicated partition 
functions is reduced to, 



FIG. 3: Several clusters for the ±J three-body Ising model on 
the triangular lattice. These are chosen to cover the whole lat- 
tice. A cluster is the unit plaquette encircled by white spins. 
The spins marked black on the original lattice are traced out 
to yield interactions among white spins in clusters. The case 
denoted by To and Uo gives the principal Boltzmann factor 
used in the naive approach. 



multicritical points in quenched random systems are thus 
different in general. In addition, the symmetry (e.g., 
spin-reversal) is quite different between two-body and 
three-body Ising models, while the above estimation is 
also given for the case of the two-body ± J Ising model 
on the self-dual lattices @, 0] • The above naive approach 
by the analysis of two face Boltzmann factors defined on 
a single triangle would not yield correct answer which re- 
flects on the detailed property of systems, the symmetry 
and the shape of the lattice under consideration. Indeed 
the detailed analysis through exact calculations by the 
virtue of the hierarchical lattice has shown deviations be- 
tween the results obtained by the above naive approach 
and the answers [fl. Ho|. It is, however, very difficult to 
obtain a precise value of the location of the multicritical 
point, that is the accuracy threshold, by numerical ap- 
proaches because of the necessity of very long time for the 
equilibration and statistical error for the average of the 
quenched randomness. In the following section, instead 
of the numerical approach, we apply an analytical theory, 
the improved method for a precise location of the multi- 
critical point [HI . The purpose of the following analysis 
is to find deviations from the p c — 0.889972 given by the 
naive approach and to detect a slight difference between 
the location of the multicritical points for the triangular 
and Union- Jack lattices. 



IV. IMPROVED METHOD 

As shown in Figs. [SJandHl let us consider to sum over 
a part of the spins, called a cluster below, on the trian- 
gular lattice and the Union- Jack lattice. The previous 
technique is only consideration of the "local" principal 
Boltzmann factors Aq and Ajj defined on an elementary 
triangle, which does not necessarily reflect the effects of 
frustration inherent in spin glasses. Then the set of 



ZM(4'UW ...) 

= zP(Af\Af\...) 



(17) 



where the superscript r stands for the type of the cluster 
that one chooses. The quantity A k is the local Boltz- 
mann factor including many-body interactions generated 
by summation over spins marked black in Fig. [3J We 
define the principal Boltzmann factors ^4q and its dual 
A*^ as those with all spins surrounding the cluster in 
the up state. We assume that a single equation gives 
the accurate location of the multicritical point for any 
number of n, similarly to the naive approach, 



A 



(r) 



A 



(18) 



This is the improved method to predict a location of the 
multicritical point with higher precision than the naive 
approach. 

One example is the hexagonal cluster as denoted by Ti 
in Fig. [3] The principal Boltzmann factors of the type 
Ti are given as 



i (TO 



n 6 

Enn^° 

{ CT «}a=lI=l 

2 n cosh™ r£2 Kt ^J 



(19) 




FIG. 4: Several clusters for the ±J three-body Ising model 
on the Union- Jack lattice. 
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Type 


# 


Pc 


I—Pc 


Slope 


T ,Uo 


1 


0.889972 


0.110028 


3.41421 


Ti 


6 


0.890212 


0.109788 


3.39884 


T 2 


18 


0.890321 


0.109679 


3.39031 


T 3 


24 0.890358 


0.109617 


3.38478 


T 4 


16 


0.890344 


0.109656 


3.38507 


T 5 


24 0.890404 


0.109596 


3.37558 


T 6 


42 


none 


none 


3.36985 


T 7 


54 


none 


none 


3.36796 


T 8 


60 


none 


none 


3.36736 


Ui 


4 


0.890725 


0.109275 


3.33658 


u 2 


16 


0.890394 


0.109606 


3.37078 


U 3 


16 


0.890825 


0.109075 


3.32365 


u 4 


36 


none 


none 


3.30885 


U 5 


64 


none 


none 


3.29998 


Perturbation 








3.20911 a 



'References Q]| and [20 



TABLE I: Location of the multicritical point of the ± J three- 
body Ising model on the triangular lattice (from Ti to T5) 
and on the Union- Jack lattice (from Ui to U3). The number 
of the triangles in the cluster is written in the second column. 
We estimate the values for the slope at the critical point T c 
and set them in the rightmost column. 



and 



.4 



*(Ti) 



n 6 

^eihi^'+^x) 

{o .a }a= l J=1 



2 3 " ^cosh 6 K + J[ri sinh 6 K^j 



> (20) 



where 77 is defined on each plaquette as shown in Fig. [3] 
We estimate a more precise location of the multicritical 
point of the ±J three-body Ising model on the triangular 
lattice as p c = 0.890212 by Eq. l[18]). 

One needs to calculate two local principal Boltzmann 
factors by summation over the internal spins, and eval- 




FIG. 5: Supplementary clusters for the computation of the 
slope of T c at p = 1. 



uate the configurational-averaged values for them as in 
Eqs. ([12]) and (|2"0"|) for the improved method. The compu- 
tational complexity of the estimation of the multicritical 
point by the improved method is 0(2 JV< )+N = '), where 
is the number of interactions in the cluster and 

(r) 

Ns is that of internal sites. In reasonable time, we per- 
form the improved method for the location of the multi- 
critical point on the triangular lattice from N^ 1 ' = 6 
to iV( T5 ) — 24, and on the Union- Jack lattice from 
AT( u = 4 to N^^ = 16 as in Figs. © and Q. The re- 
sults obtained by the improved method are summarized 
in TableQ] As you can see, all the results by the improved 
method show higher values p c than p c — 0.889972 by the 
naive approach. In the context of the quantum error cor- 
recting code, the error threshold of the topological color 
code is lower than the probability 1 — p c = 0.110028 for 
the quantum Gilbert- Varshamov bound with zero rate 
encoding as in 

Table EllH Eft 

Next let us to examine the performance of the im- 
proved method, how close to the exact answer, for the 
random three-body Ising model on the triangular and 
Union- Jack lattices. If we formulate the above analysis 
without restriction to the Nishimori line, we can give an 
approximative shape of the phase boundary [Iol | . We es- 
timate the value of the slope at the critical point T c for 
p = 1 as listed in Table [IJ It is one of the references 
for the performance of the improved method to compare 
these values with the solution obtained by a simple per- 
turbation 3.20911 [lj|, which is essentially the same way 
as Domany's result for the two-body Ising model on the 
square lattice [2(|. Since the computational complexity 
can be reduced to 0(N^2 N = ' ) for the calculation of the 
slope of T c at p = 1 by the improved method, we eval- 
uate only the values of the slope by using larger size of 
the clusters denoted by Tg, T7, and T§ as in Fig. [5] As 
you can see in Table. [T] for the case for the triangular 
lattice, if the number of the interactions in the cluster 
becomes larger, the estimations for the slope of T c by the 
improved method steadily decrease into the perturbation 
solution 3.20911 but the convergence into the value is ex- 
tremely slow. Such a behavior of the slope of T c at p = 1 
is a sharp difference from the case for the two-body ± J 
Ising model on the square lattice [To| . Though it is thus 
expected that the predictions of the location of the multi- 
critical point by the improved method certainly approach 
the exact answer, our estimation p c — 0.8902 — 4 for the 
random three-body Ising model on the triangular lattice 
is a delicate conclusion. At least, we can state that the 
accuracy threshold 1— p c for the color code on the hexago- 
nal lattice is lower than the quantum Gilbert- Varshamov 
bound with a zero rate encoding. 

On the other hand, we obtain a relatively reliable result 
for the Union-Jack lattice, though we find an exception of 
the systematic improvement by the use of larger size clus- 
ters the result by U2 cluster. From the difference between 
two values for the slope of T c at p — 1 by the improved 
method and the perturbation solution, U2 cluster is con- 
sidered to be an wrong approximation for estimations of 
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the location of the multicritical point. It is considered 
that the shape of U2 cluster is slightly different from the 
series of other clusters as shown in Fig. |4j If we use a 
series of the similar clusters to Ui and U3 (see Fig. [5]) for 
the further approximations for the slope of T c at p = 1 , 
we find convergence into the perturbation result 3.20911, 
which can be derived in the same way in Refs. and 
[20T |. with a similar degree to the case for the two-body 
±J Ising model on the square lattice [l0| . We thus pre- 
dict that the location of the multicritical point for the 
Union- Jack lattice would be p c — 0.8907 — 8 from the 
results by Ui and U3 clusters, or equivalently the accu- 
racy threshold for the color code on the square-octagonal 
lattice is 1 - p c = 0.1092 - 3. 

The improved method thus strongly depends on the 
shape of the used cluster, as well as the number of in- 
teractions included in the cluster. It is considered that 
the systematic improvement with convergence into the 
answer is yielded by the ingenious choice of the shape of 
the cluster following properties and symmetry of the lat- 
tice under consideration. We would be able to approach 
the exact answer when we choose the clusters keeping the 
underlying symmetry in mind in a systematic way. The 
slow convergence for the case on the triangular lattice 
would be caused by an inadequate choice of the clusters 
on it. This point should be clear in the future. 

V. CONCLUSION 

We predicted the accuracy thresholds of the topolog- 
ical color codes on the hexagonal and square-octagonal 
lattices, which are kinds of the quantum error correcting 
codes exploiting the topological properties in systems. 
These accuracy thresholds correspond to the locations of 
the multicritical points for the random three-body Ising 



models on the triangular and Union-Jack lattices, respec- 
tively. The accuracy thresholds of these two topological 
color codes are expected to be lower than the probabil- 
ity for the quantum Gilbert- Varshamov bound with zero 
rate encoding. Accepting our estimations, we find that 
the accuracy thresholds of two color codes are equiva- 
lent to that of the quantum toric code on the square 
lattice, which corresponds to the location of the multi- 
critical point for the ± J Ising model on the square lattice 
[l|. The accuracy threshold is given as 1— p c = 0.1092 — 3 
from several results by the improved method 1— Pc = 
0.10919(7) by highly precise Monte-Carlo simulation [2l| . 
and a very recent estimation 1 — p c = 0.10939(6) by the 
concept of the conformal field theory (2^]. In conclu- 
sion, the topological color code on the hexagonal and 
square-octagonal lattice would have equivalent robust- 
ness to that of the topological toric code on the square 
lattice despite their computational capabilities. This fact 
is indeed shown by the numerical simulation performed in 
very recent work by Katzgraber et al. @. The present 
work provided the same conclusion but would be valu- 
able as the analytical predictions of the topological color 
codes on the hexagonal and square-octagonal lattices. 
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